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PACS 72. 10. -d - Theory of electronic transport; scattering mechanisms 
PACS 73 . 20 . At - Surface states, band structure, electron density of states 
PACS 75 . 70 . T j - Spin-orbit effects 

Abstract - We investigate the effects of spin-flip scattering on the Hall transport and the spectral 
properties of gapped Dirac fermions. We find that in the weak scattering regime, the Berry 
curvature distribution is dramatically compressed in the electronic energy spectrum, becoming 
singular at band edges. As a result the Hall conductivity has a sudden jump (or drop) of e 2 /(2/i) 
when the Fermi energy sweeps across the band edges, and otherwise is a constant quantized in 
units of e 2 /(2h). The spectral properties such as the density of states and the spin polarization 
are also greatly enhanced near band edges. Possible experimental methods to detect these effects 
are discussed. 
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Introduction. — Dirac fermion has found its ubiq- 
uitous appearance in interesting physical systems such as 
graphene [JJ and topological insulators [2)15] . The most im- 
portant feature of Dirac fermion is its strong (pscudo)spin- 
orbit coupling, which is directly encoded in the Hamilto- 
nian and underlies many unusual effects. When certain 
symmetry breaking mechanism is introduced, a gap can 
be opened up at the Dirac point. Gapped Dirac fermion 
(GDF) is in some sense more interesting in that its elec- 
tronic band develops a Berry curvature gauge field [HEI 
which acts like a magnetic field in the momentum space, 
and deflects electron flow in transverse directions leading 
to the anomalous Hall effect [BJ. 

Disorder scattering usually has strong effects on the Hall 
transport even in the weak disorder limit when the system 
is in a metallic state. So far, most studies of scattering 
effects on Dirac fermions have focused on spin independent 
scattering. However, due to the strong spin-orbit coupling, 
it is natural to expect that the carrier motion should have 
a sensitive dependence on the change of spin state during 
a scattering. 

In this paper, we study the effects of spin-flip scatter- 
ing on the properties of GDFs. Our work is motivated by 
recent advances in creating GDF on the surface of topo- 
logical insulators 08] and the ability of systematic control 
of surface carrier density through doping or gating [5] . We 
show that in the weak scattering regime spin-flip scatter- 
ing compresses the Berry curvature distribution dramati- 
cally towards the band edge. Hence whenever the Fermi 



energy is tuned across a band edge, the Hall conductivity 
has a sudden change of e 2 /(2h), i.e. half of the conduc- 
tance quantum. Away from the band edges, the Hall con- 
ductivity is a constant independent of the carrier density. 
Since Hall transport is central to several exotic physical ef- 
fects [10] proposed recently for topological insulators, the 
discovery presented here is expected to have important 
observable consequences. We also find that the spectral 
properties like the density of states (DOS) and the spin 
polarization which are usually not very sensitive to dis- 
order scattering get greatly enhanced near band edges by 
the spin-flip scattering, which can be detected by photoe- 
mission spectroscopy or tunneling spectroscopy measure- 
ments. 

Quantized Hall conductivity and Berry curva- 
ture compression. — A GDF is described by the Hamil- 
tonian 

Ho = V F k -T + At z , (1) 

where vf is the Fermi velocity, k = (k x ,k y ,Q), and 
r = (t x ,t v ,t z ) is a vector of Pauli matrices acting on 
the spin degrees of freedom. Without the second term, 
the Hamiltonian is for a massless Dirac fermion with a 
cone-like energy dispersion. The term At 2 induces a finite 
mass and opens up an energy gap of 2A in the spectrum. 
Such a term can be generated by breaking the time re- 
versal symmetry at the surface of a topological insulator, 
for example, through magnetic doping [7] or coating with 
insulating ferromagnetic films. 
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Because of the strong spin-orbit coupling and the spin 
splitting induced by the mass term, an anomalous Hall ef- 
fect in the absence of external magnetic field is expected. 
In the weak scattering regime, the anomalous Hall con- 
ductivity has an important contribution cr® y from the mo- 
mentum space Berry curvature of the spin-orbit coupled 
bands [BJ. It is known as the intrinsic contribution be- 
cause it is an intrinsic property of a crystal. For a two- 
dimensional (2D) system the intrinsic contribution is given 
by the integral of momentum space Berry curvature fl n k 
of all the occupied states \nk) (n is the band index), 

a % = -iiE„)e ttnkfnk, where A is the area of the 2D 
system, is the Fermi distribution function, 

^ 2Im(nk\v x \n'k)(n'k\vy\nk) , n . 

Vnk = - > , — \o~ , (2) 



(a) 



Vi (= VFTi/h with i = x,y for GDF) is the velocity oper- 
ator and hu} n k is the energy of the state \nk). The dis- 
tribution of Berry curvature in the energy spectrum can 
be described by the density of Berry curvature 51(e) = 
I Enfe ^nit^ — Fv^nk)- Then the intrinsic contribution 
can be put into the form a xy = — ^- J ded(e)f(e). 

For GDFs, the density of Berry curvature and the in- 
trinsic contribution of Hall conductivity can be easily cal- 
culated using the formulae above. The result is 
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sgn( e )6(| £ | - |A|), 
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e(M-|A|)+ B gn(A)e(|A|-M 



(4) 

where Q(x) is the Hcaviside step function and Sf is the 
Fermi energy. As shown by the blue curves in fig. [TJ while 
the density of Berry curvature has its maximum (or min- 
imum) value at band edges, its distribution spreads over 
the whole energy spectrum, hence the value of a xy changes 
gradually as a function of the Fermi energy. For each band 
the total weight of Berry curvature corresponds to a con- 
tribution of ±e 2 / (2/i) to the Hall conductivity, where the 
sign difference reflects the different helicities of the two 
bands. It needs to be mentioned that the half quantized 
contribution from the filled lower band does not contra- 
dict the well known fact that the contribution from a com- 
pletely filled band must be integer quantized [BJ. This is 
because the Dirac model is energetically unbounded. It 
can only serve as a low energy effective model for any real 
physical systems, where the evaluation of the contribution 
from all the completely filled bands obviously has to go be- 
yond the effective model. On the other hand, it is crucial 
to notice that the disorder related contribution, which is 
our focus in this paper, is tied to the Fermi surface hence 
is well captured within such effective model. 

The existence of disorder scattering is essential for sta- 
bilizing the electron distribution under a drive field and 
establishing the dynamic steady state. In the case of the 
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Fig. 1: (color online), (a) The anomalous Hall conductivity 
and (b) the density of Berry curvature for the gapped Dirac 
fermion (with A > 0) as a function of energy. The blue curve 
in (a) is for the intrinsic contribution and the red curve is for 
a disordered system with weak spin-flip scattering. The Berry 
curvature is compressed by the scattering to be singular at 
band edges, as shown schematically in (b). 



anomalous Hall effect, the role of disorder scattering is es- 
pecially important yet complicated. In the weak scattering 
regime, when expanded in terms of the disorder density 
ridis) the leading contribution to cr xy is generally of or- 
der n7 ; , known as the skew scattering contribution which 
originates from the asymmetry in the scattering rate and 
depends on the details of the disorder model [BJ . The next 
leading order (of n^ ls dependence) consists of the intrinsic 
contribution discussed above and the side jump contribu- 
tion [BJ. The skew scattering contribution and the side 
jump contribution arise from disorder scattering, they are 
referred to as extrinsic contributions. The side jump con- 
tribution is peculiar in that it arises from scattering but is 
independent of both the scattering strength and the disor- 
der density. Due to its presence, even for disorder models 
with vanishing skew scattering (such as Gaussian disorder 
models), the Hall conductivity a xy in the weak disorder 
limit rcdis — * + docs not correspond to the value of the 
intrinsic contribution a xy [BJ. This is an important point 
that needs to be emphasized. 

In a recent study of the anomalous Hall effect [11] , we 
find that the extrinsic contributions to a xy depend sensi- 
tively on the spin structure of disorder scattering. This 
is understandable because the Hall transport is a result 
of spin-orbit coupling, different actions on the spin state 
during a scattering will also strongly affect the carrier's 
orbital motion. Here we are most interested in the ef- 
fects of spin-flip scattering on the gapped Dirac fermions. 
Such kind of scattering can be modeled as a random po- 
tential U(r) = J2 a Ba ( r -Ra)-r = B x (r)T x + B y (r)r y , 
with B a being a random in-plane vector. We assume that 
the random field satisfies the correlation (Bi(q)B*(q')) c = 
n<ii s u 2 5ij5(q — q'), where i, j £ {x,y}, Bi(q) is the Fourier 
component of Bi(r), (■ ■ -) c stands for disorder average and 
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u is the scattering strength. We assume that the system 
has no preferred in-plane direction, such that the random 
vector B a is uniformly distributed in plane. Then the 
third order disorder correlation (UUU) C must vanish iden- 
tically. Because the skew scattering contribution depends 
on the third order correlation [6]|T2, this means that the 
skew scattering process is forbidden for this kind of disor- 
der. Therefore the leading order contribution to the Hall 
conductivity comes from the intrinsic and the side jump 
terms which are independent of scattering strength and 
disorder density. 

In the weak scattering regime, the Hall conductivity can 
be evaluated perturbatively using the Kubo-Streda for- 
mula Q2] ■ The calculation has been detailed in our previ- 
ous work [11] . Here we quote the final result which is that 
the side jump contribution is 
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hence the total Hall conductivity is given by 
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This is a very interesting result. As plotted in fig. [lja), 
a xy has a sudden jump or drop of half the conductance 
quantum at the band edges. As soon as the Fermi energy 
passes the band edge, the value of a xy becomes a constant 
quantized in units of e 2 /2h. The effect of the spin-flip 
scattering makes the anomalous Hall conductivity quan- 
tized in the whole energy spectrum, which is different from 
other types of disorders [TTHT2] . 

In this paper, we would like to promote another per- 
spective on the unusual phenomena described by eq. (pj). 
Instead of starting from a clean system and treating disor- 
ders as perturbations, we start directly from a disordered 
system and try to generalize the concept of Berry curva- 
ture for the disordered case. For each cigenstate \a) with 
eigen-energy huj a of the disordered system, we define its 
Berry curvature as 
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We further introduce a Berry curvature operator SI = 
J2 a fla\<x){&\, such that the density of Berry curvature 
can be generalized as 



(8) 



It follows from the Kubo formula that the Hall conductiv- 
ity for disordered systems can be expressed as 



-- / de(n(e)) c f(e), 



(9) 



with a disorder averaged density of Berry curvature. This 
is an exact result. We again emphasize that although 
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Fig. 2: (color online). Numerical result of (a) the Hall con- 
ductivity and (b) the Berry curvature distribution for in- 
creasing system sizes with Nk = 40, 50, 60, 70 averaged over 
200, 100, 80, 50 samples respectively. In the calculation we 
set vf = A = 1, the momentum cutoff K c = 6, and the 
strength |B,(q)| (i — x,y) is a random number uniformly 
distributed in an interval [0, W] with W = 0.03 (accordingly 
ridisU 2 ~ W 2 K 2 /3 ~ 0.01). The error bar stands for one stan- 
dard deviation. 



the definition of generalized Berry curvature recovers that 
for a perfect crystal when the eigenstatcs \a) are simply 
Bloch states \nk), the weak disorder limit (ridis — > + ) of 
eq. © is not equal to the intrinsic contribution. It also 
contains the side jump contribution (for vanishing skew 
scattering). Mathematically, this is because the disorder 
averaged Berry curvature contains nontrivial vertex cor- 
rections which do not vanish in the weak disorder limit. 
From eq. ((5]) we find that in the weak disorder limit, 
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As shown schematically in fig. [Tfb) , the density of Berry 
curvature, which is originally spread over the energy spec- 
trum, gets compressed dramatically towards the band 
edges by the scattering, while its total weight (±^r for 
lower and upper band respectively) is kept unchanged. 

Based on eqs. d)-© we perform a numerical calcula- 
tion to further confirm our analytic result. Due to the 
particle-hole symmetry, a xy is a symmetric function with 
respect to Ef = 0, hence here we only focus on the positive 
energy range. The calculation is performed in momentum 
space on disordered systems of size Nk x Nk ■ Figure [5] 
shows the results for a xy and SI with different system size 
Nu ■ The physically meaningful value of a transport coeffi- 
cient such as a xy corresponds to the thermodynamic limit 
when Nk — > oo. The result in fig. [5] indeed shows that the 
density of Berry curvature is squeezed by the weak spin- 
flip scattering to the upper band bottom. Consequently 
the contribution to the Hall conductivity from the upper 
band approaches a half-quantized plateau in the energy 
spectrum, which confirms our analytic result. 

To completely understand the microscopic mechanism 
of this peculiar disorder effect is difficult. Nevertheless, 
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we notice that for the gapped Dirac fermion, the spin-flip 
scattering has weak scattering rates near the band bot- 
tom because states there have large spin z-components, 
meanwhile large scattering rates occur for high energies 
where the spins are more or less lying in the plane. Our 
result indicates that the Berry curvature tends to be ex- 
pelled from the strong scattering region towards the weak 
scattering region while keeping its total weight unchanged, 
hence making the curvature distribution squeezed towards 
band bottom. This is similar to the migration of Chcrn 
number carrying states driven by impurity scattering ob- 
served in the quantum Hall effect [T4"l[l"5] . 

Enhanced DOS and spin polarization. Spectral 
properties such as the DOS are usually considered to be 
insensitive to weak disorder scattering. This is particularly 
true when the DOS does not vary much in the energy 
range under consideration, such as for a 2D free electron 
gas where the DOS is a constant. A GDF gas without 
disorder has a DOS 
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There are three salient features of this DOS that are cru- 
cial for the effects to be discussed below: (1) p increases 
linearly from the band edge; (2) p has a finite value ij^r 

at band edges; (3) po does not depend on the value of the 
spin splitting A as long as |e| > |A|. Figure [SJa) shows the 
numerical result of the DOS for a clean system together 
with the DOS for a disordered system with spin-flip dis- 
orders. It is observed that spin- flip scattering strongly 
enhances the DOS near the band edge, and slightly de- 
creases the DOS at large energies, as being required by 
the conservation of the total number of states. 

This effect can be qualitatively understood from the 
self-energy correction due to scattering. In the Born 
approximation, (assuming e > A > in the follow- 
ing), the 2x2 retarded self energy is given by Y> R = 
^E k (UG§(k,e)U) Cl with 
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and Im£ fl (e) = — (er - At z ), where e c is the en- 
ergy cutoff and r is the 2x2 identity matrix. While 
the imaginary part of the self-energy only contributes to a 
small level broadening in the weak scattering regime, the 
real part has a large magnitude near the band edge and 
diverges logarithmically as e — > A. In fact, this disorder 
induced singularity is a generic feature for systems having 
a finite (unperturbed) DOS at the band edge, as is the 
case for the GDF system. 

Thermodynamic properties can be extracted from the 
Green function G R = (e-Hq-^)- 1 . From cq. $FZ$, the 
self-energy only corrects the energy argument e and spin 
splitting A in G R . It follows that the disordered DOS can 
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Fig. 3: (color online). Numerical results of (a) the DOS and 
(b) the spin polarization S z for increasing disorder strength 
with W = 0.02, 0.04, 0.06 with each data point averaged over 
50 samples, where W is defined in the caption of fig. [2] In the 
calculation, we set Nt = 70, and Vf = A = 1. 



be approximated as 

p( £ )~ /jo[e(l + A),A(l + A)], (13) 

with the right hand side standing for the DOS at en- 
ergy e(l + A) for a clean GDF system with A replaced 

by A(l + A), where A = In ^ is the prefactor 



appearing in eq. (|T2"T) . Near the band edge, A is a posi- 
tive number, hence p(e) ~ (1 + A)po(s) is enhanced by a 
factor of (1 + A) compared with the clean case [IB]. Now 
it is clear that a strong enhancement of DOS by scatter- 
ing requires the original DOS has a large variation on the 
energy scale of A A. Therefore such effect does not occur 
for a conventional 2D free electron gas, and although the 
GDF near the band edge resembles the 2D free electron 
gas, the self-energy correction makes the Dirac behavior 
at higher energies relevant. Similarly, the spin polariza- 
tion S z (e) = —^jlmTi[T z G R (e)] is also enhanced at band 
edge by the factor (1 + A), as shown in fig. Elb). 

Discussion. — The exotic effects reported here arise 
from the interplay between the GDF and the spin-flip scat- 
tering. There are two important properties of GDFs un- 
derlying these effects: the Berry curvature and the un- 
usual DOS. Disorder scattering modifies these fundamen- 
tal properties in the spectrum, thereby strongly influences 
the system behavior. 

Spin conserving scattering such as B (r)T or B z (r)r z 
also affects the anomalous Hall transport of GDF system 
in metallic state. In contrast with the spin-flip scatter- 
ing, their skew scattering contribution is in general non- 
vanishing. Their side jump contributions are also distinct 
from that of the spin-flip scattering [Tl][T2]. In particu- 
lar, the Berry curvature is not compressed into a small 
region in the spectrum. The combined effects of different 
kinds of scattering would depend on the competition be- 
tween them fTT1ll7j . As for the disorder modified DOS, the 
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spin conserving scattering has similar effect as the spin-flip 
scattering. This can be understood by noticing that the 
self-energy correction for spin conserving scattering only 
differs from eq. (fT2|) by a sign change of A. Because po 
does not depend on A (for energies within the bands) , the 
result is similar to that for the spin-flip scattering. Mean- 
while, the spin polarization becomes S z (e) ~ (1 — X)S§(e). 
In energy range where A > 1, the spin polarization can 
even have a sign change. 

GDF has been realized on the surface of topological in- 
sulators [7|. This is partly motivated by the possibility 
to achieve a half-quantized anomalous Hall effect by tun- 
ing the Fermi level into the surface band gap [51fT51. In 
contrast, our focus here is instead on the energy range out- 
side the gap. We notice that for a purely 2D system with 
embedded gapped Dirac spectrum, it is possible to real- 
ize half-quantized anomalous Hall effect in the metallic 
state with the help of spin-flip scattering. For example, 
consider a Hamiltonian Hq = vp (sin k x r x + smk y T y ) + 
(A + 2 — cosk x — cosky)T z (A > 0) which corresponds to 
a tight-binding lattice model proposed to model the Mn 
doped HgTe/CdTe quantum wells [TDJ. It has a direct gap 
at F point of the Brillouin zone. The low energy effec- 
tive model near the gap is just that for a GDF. As we 
mentioned before, the contribution to the Hall conductiv- 
ity from the completely filled lower band (or any remote 
valence band) must be integer quantized. Therefore once 
the Fermi energy is slightly above the upper band bottom 
and spin-flip scattering is turned on, an e 2 /(2h) contribu- 
tion arises from the upper band bottom making the total 
Hall conductivity half-quantized. In fact, the non- integer 
value of Hall conductivity dictates that the system must 
be conducting [2D]. The effect itself would be very inter- 
esting since it realizes a fractional quantized Hall effect for 
non-interacting electrons in a metallic state, plus it does 
not require the completely filled valence bands to have a 
nontrivial topology, i.e. with a nonzero Chern number [5]. 

Finally, we briefly comment on the possible experimen- 
tal methods to detect the effects studied here on the sur- 
face of topological insulators. The surface band gap can 
be opened by depositing an ultrathin insulating ferromag- 
netic film on top of a topological insulator such that a 
proximity-induced exchange coupling can be introduced. 
Spin-flip scattering can be realized by thermally excited 
spin waves in the ferromagnetic film. However, to make 
spin wave scattering dominate over other scattering pro- 
cesses is a nontrivial task. It requires a temperature range 
where spin wave excitation dominates over phonon exci- 
tation and improved sample quality such that very few 
impurities are present near the interface. Fcrromagetic 
insulating materials such as EuS and EuO with the Curie 
temperature much lower than the Debye temperature are 
possible candidates. The Hall conductivity can be mea- 
sured through standard multi-terminal setup. The sudden 
change of Hall conductivity at band edges will also man- 
ifest in the signal of magneto-electric or magneto-optic 
effects [TDJ . Enhanced DOS and spin polarization at band 



edges can be directly detected through the photoemission 
spectroscopy or the tunneling spectroscopy measurements, 
which should be much easier. 

* * * 

The authors thank N. A. Sinitsyn and J. Sinova 
for helpful discussions. This work is supported by 
NSFC (10974231) (Y.Y.), DOE (DE-FG03-02ER45958) 
and Texas Advanced Research Program (Q.N.). 

REFERENCES 

[I] Castro Neto A. H. et al, Rev. Mod. Phys., 81 (2009) 
109. 

[2] Fu L., Kane C. L. and Mele E. J., Phys. Rev. Lett, 98 
(2007) 106803; More J. E. and Balents L., Phys. Rev. 
B, 75 (2007) 121306(R); Qi X.-L., Hughes T. L. and 
Zhang S.-C, Phys. Rev. B, 78 (2008) 195424; Hsieh D. 
et al, Nature, 452 (2008) 970; XlA Y. et al, Nat. Phys., 
5 (2009) 398; Chen Y. L. et al, Science, 325 (2009) 178. 

[3] Hazan M. Z. and Kane C. L., Rev. Mod. Phys., 82 (2010) 
3045. 

[4] Bohm A., Mostafazadeh A., Koizumi H., Niu Q. and 
Zwanziger J., The Geometric Phase in Quantum Systems 
( Springer- Verlag, Berlin) 2003, p. 19. 

[5] Xiao D., Chang M.-C. and Niu Q., Rev. Mod. Phys., 82 

(2010) 1959. 

[6] NAGAOSA N. et al, Rev. Mod. Phys., 82 (2010) 1539. 
[7] Chen Y. L. et al, Science, 329 (2010) 659. 
[8] Wray L. A. et aL larXiv:1009.6216l 
[9] XlA Y. et al.. larXiv:0907.3089l 

[10] Garate I. and Franz M., Phys. Rev. Lett., 104 (2010) 
146802; Tse W.-K. and MacDonald A. H., Phys. Rev. 
Lett, 105 (2010) 057401; Yokoyama T. et al, Phys. Rev. 
B, 81 (2010) 241410(R). 

[II] Yang S. A., Pan H., Yao Y. and Niu Q., Phys. Rev. B, 
83 (2011) 125122. 

[12] Sinitsyn N. A. et al, Phys. Rev. B, 75 (2007) 045315. 
[13] Streda P., J. Phys. C, 15 (1982) L717. 
[14] Liu D. Z., Xie X. C. and Niu Q., Phys. Rev. Lett, 76 
(1996) 975. 

[15] Sheng D. N., Weng Z. Y. and Wen X. G, Phys. Rev. 

B, 64 (2001) 165317. 
[16] A more careful evaluation leads to the result that p(e) = 

Po(s) [l + n £ a v f ( m | fiZ^-i | ~ 2 )] > whicn shows that the 

curves of p(e) and po(e) have a crossing point at the energy 
2 — 2 I 

e* when In grz^? — 2 = 0. With the parameter values in 
our numerical calculation, we find that e» = 2.55 which 
agrees very well with the result in Fig[3Ja). 
[17] Nomura K. and NAGAOSA N., Phys. Rev. Lett., 106 

(2011) 166802. 

[18] Chu R. et a UarXiv: 1007.04971 

[19] Liu C.-X. et al, Phys. Rev. Lett, 101 (2008) 146802. 
[20] Niu Q., Thouless D. J. and Wu Y.-S., Phys. Rev. B, 
31 (1985) 3372. 



p-5 



